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First we have assumed the most general static spherically symmetric black hole metric. The 
accretion of any general kind of fluid flow around the black hole have been investigated. The 
accretion of fluid flow around the modified Hayward black hole have been analyzed and we then 
calculated the critical point, fluid 4 velocity and velocity of sound during accretion process. Also the 
nature of the dynamical mass of black hole during accretion of fluid flow and taking into consideration 
of Hawking radiation from black hole i.e., evaporation of black hole have been analyzed. 
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I. INTRODUCTION 

At present we live in a Universe which is expanding 
and the expansion rate is increasing i.e, the Uni¬ 
verse is accelerating which was confirmed by recent 
Supernova type la observations 0 i- The large 
scale structure 00, cosmic microwave background 
radiation 0, WMAP observations 0-0 also support 
this acceleration of the Universe. This acceleration 
is caused by some unknown matter which produce 
strong sufficient negative pressure (with positive 
energy density), known as dark energy. The present 
Universe occupies ~ 4% ordinary matter, ~ 74% 
dark energy and ~ 22% dark matter. Since dark 
energy and dark matter are main two components 
in our universe, such that the present dark energy 
and dark matter densities are 7.01 x 10 ~ 27 kg/m 3 
and 2.18 x 10 ~ 27 kg/m 3 , respectively. Most simplest 
candidate of dark energy is the cosmological constant 
A which obeys the equation of state EoS p = wp with 
EoS parameter w = — 1 00- Another candidates of 
dark energy are quintessence (where EoS parameter 
satisfies — 1 < w < —1/3) fill . Il2| and phantom 
(where EoS parameter satisfies w < — 1) 0]. Till 
now there are lot of dark energy models have been 
considered. A brief review of dark energy models is 
found in the ref. 0- 

A condensed object (e.g. neutron stars, black 
holes, etc.) surrounded by a fluid can capture 
particles of fluid that pass within a certain distance 
from the condensed object. This phenomena is 
termed as accretion of fluid by condensed objects. 
In Newtonian theory of gravity, the problem of 
accretion of matter onto the compact object was first 
formulated by Bondi (l5j| . Michel [l6j first obtained 
analytic relativistic accretion (of gas) solution onto 
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the static Schwarzschild black hole. Such accretion 
processes are candidates to mechanisms of formation 
of supermassive black holes (SMBH) in the center 
of most active galaxies 0- In particular, it should 
follow some analogies with the process proposed by 
Salpeter et al 0] where galaxies and quasars could 
get some of their energy from processes of accretion. 
Using this accretion procedure, Babichev et al 00 
formulated the accretion of phantom dark energy 
onto a static Schwarzschild black hole and shown that 
static Schwarzschild black hole mass will gradually 
decrease due to strong negative pressure of phantom 
energy and finally all the masses tend to zero near the 
big rip singularity. Sun [0 discussed phantom energy 
accretion onto a Black hole in the cyclic universe. 
Jamil 0 has investigated accretion of phantom like 
modified variable Chaplygin gas onto Schwarzschild 
black hole. Phantom energy accretion by stringly 
charged black hole has been discussed by Sharif et 
al [23l |. Dark matter and dark energy accretion onto 
static black hole has been discussed by Kim et al 
[ 24 } . Also the accretion of dark energy onto the 
more general Kerr-Newman black hole was studied 
by Madrid et al [25). The new variable modified 
Chaplygin gas and generalized cosmic Chaplygin gas 
dark energy accretions and onto Kerr-Newman black 
hole and their features were studied Bhadra et al 0] . 
Several authors have discussed the accretions 

of various components of dark energy onto several 
types of black holes. 

In the present work, first we assume the most gen¬ 
eral static spherically symmetric black hole metric in 
section II. The accretion of any general kind of fluid 
flow around the black hole will be investigated. The 
accretion of fluid flow around the modified Hayward 
black hole will be analyzed in section III and we then 
calculate the critical point, fluid 4 velocity and ve¬ 
locity of sound during accretion process. Also the 
nature of the dynamical mass of black hole during 
accretion of fluid flow and taking into consideration 
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of Hawking radiation from black hole i.e., evaporation 
of black hole will be analyzed in section IV. Finally 
we shall draw fruitful discussions about accretion of 
fluids upon modified Hayward black hole in section V. 


II. ACCRETION PHENOMENA OF 
GENERAL STATIC SPHERICALLY 
SYMMETRIC BLACK HOLE 

First we consider general static spherically symmet¬ 
ric metric given by 

ds 2 = — A(r)dt 2 H-— dr 2 + r 2 (d0 2 + sinOdq i> 2 ) (1) 

B(r) 

where A(r) > 0 and B(r ) > 0 are functions of r only. 
We can choose A(r) and B(r) in such a way that the 
above metric represents a black hole metric. Let us 
assume M is the mass of the black hole. For instance, 
if A(r) = B(r ) = 1— ; the above metric represents 

Schwarzschild black hole. 

The energy-momentum tensor for the fluid is given 

by 


= (p + p)up.u v + pg^ 


( 2 ) 


where p and p are the energy density and pres¬ 
sure of the fluid. The four velocity vector of the 


fluid flow is given by u p = 


dx* 


= (w, it 1 ,0,0) 


where u° and u 1 are the non-zero components of 
velocity vector satisfying = —1. This im¬ 

plies goou°u° 

1,2 . 


gnu 1 u 1 


= — 1. So we can obtain 


( u °) 2 = \b B an b let the radial velocity of the flow 

u 1 = u , so we have uq = goou 0 = ' A - 


B. Here 


y/—g = r 2 sin6. From above equation ©, we 

obtain Tq = (p + p)uqu. It is assumed that u < 0 for 
inward flow of the fluid towards the black hole. 


flux equation can be derived by the projection of the 
conservation law for energy-momentum tensor onto 
the fluid four-velocity, i.e., UpTfp' = 0, which gives 
u^p ifJj + (p + p)uft = 0. From this, we obtain 


ur 2 M~ 2 


[A 

r r ph dp 

\l— exp 

.Jp a c p+p(p). 


= —C 


(4) 


where C is integration constant (energy flux onto the 
black hole) and the associated minus sign is taken for 
convenience. Also ph and p^ represent the energy 
densities at the black hole horizon and at infinity re¬ 
spectively. Combining equations © and ©, we ob¬ 
tain, 


{p + p)y/vF' 


[A 

r ph dp 

\l — exp 

. Jpo. P+P(p). 


= C 2 (5) 


where, C 2 = —C\/C = p^ + p{poo)■ The equation of 
mass flux Jft = 0 is given by ^ = 0’ which 

integrates to pu l y/—g = A\ and yields 


pur 2 M 2 



( 6 ) 


where, C3 is the integration constant. From © and 
©, we obtain, 


p ~b p A / ~ ~ C\ „ 

- \ — y u 2 + B = — = Ci= constant 


B 


C* 


Now let us assume, 


2 = dhi(p + p) _ 
din p 

So from equations ©, © and ©, we obtain 


V 2 - 


B 


du 

u 


-2V 2 


1 (A 


A 


B' 

~B 


(7) 

( 8 ) 


In the fluid flow, we may assume that the fluid is 
dark matter or any kind of dark energy. A proper 
dark-energy accretion model for static spherically 
symmetric black hole should be obtained by gener¬ 
alizing the Michel’s theory |l6j. In the dark energy 
accretion onto Schwarzschild black hole, Babichev et 
al [ljj, [20| have performed the above generalization. 
We shall follow now the above procedure in the case 
of static spherically symmetric black hole. The rela¬ 
tivistic Bernoulli’s equation (the time component) of 
the energy-momentum conservation law T.^ u = 0, we 
obtain ^ (Tg 1 yj—g) = 0 which provides the first in¬ 
tegral, (p + p)uou 1 y/^g = Ci, that simplifies to 

ur 2 M~ 2 (p + p)^\/u 2 + B = Ci (3) 

JD 


x(P 2 + l )r + 


rB' 

2 (u 2 + B) 


*=0 

r 


(9) 


Now if one or the other of the bracketed terms in © 
vanishes, we get a turn-around point and in this case, 
the solutions will be the double-valued in either r or u. 
There are only solutions which pass through a critical 
point that correspond to material falling into (or flow¬ 
ing out of) the object with monotonically increasing 
velocity along with the particle trajectory. A point 
where speed of flow is equal to the speed of sound, 
such a point is called a critical point. It is assumed 
that the critical point of accretion is located at r = r c 
which is obtained by taking the both bracketed terms 
(coefficients of du and dr) in Eq. © to be zero. So 
at the critical point, we obtain 


where Ci is the integration constant which has the 
dimension of the energy density. Moreover, the energy 


W 2 = 


ul + B(r c ) 


( 10 ) 
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and 

4K 2 


A'(r c ) B'(r c ) 


A(r c ) B(r c ) _ 


(V 2 + 1) + B '^ (11) 

1 c + j + u 2 + B(r c j> 


Here, subscript c denotes the critical value and u c is 
the critical speed of flow at the critical point r c . From 
above two expressions, we have 


u 2 = 


B'(r c ) A'(r c ) 


and 
V 2 = 


1 + 2 


2 A(r c ) 

A(r c ) B(r ) 


A'(r c ) B'(r c ) 
A(r c ) B(r c ) 


1 -1 


( 12 ) 


A\r c ) B'(r c ) 
A'(r c ) B'{t c ) \r c A(r c ) B(r c ) 


(13) 


At the critical point r c , the sound speed can be de¬ 
termined by 


2 _ d P 
° s dp 


C A V c (V 2 + 1) B(r c ) 

u c V A(r c ) 


-1 (14) 


The physically acceptable solutions of the above equa¬ 
tions may be obtained if u 2 > 0 and V 2 > 0 which 
leads to 

A'(r c )B'(r c ) > 0 and 1 > 4^ - (15) 


From the above equation we can obtain the bound of 
r c if A and B are known for several kinds of static 
black holes. 


III. ACCRETION PHENOMENA OF 
MODIFIED HAYWARD BLACK HOLE 

The static spherically symmetric space-time is de¬ 
scribed by the Hayward metric which is obtained by 
A(r) = B(r) in equation ([1|) and is given by [35j 

ds 2 = —B(r)dt 2 H- - dr 2 + r 2 (d0 2 + sinddep 2 Y16) 

B(r) 

Here, M is the mass of Hayward black hole and 
B(r) = 1 — r 3+!> r Ml 2 i where l is a parameter with 
dimensions of length (Hubble length) with small scale 
related to the inverse cosmological constant A (l is 
a convenient encoding of the central energy density 
gAp ~ A, assumed positive). Such behavior has been 
proposed by Sakharov [36!, L371 as the equation of state 
of matter at high density and by Markov [38|, [39} based 
on an upper limit on density or curvature, to be ul¬ 
timately justified by a quantum theory of gravity. In 
the limit r —> oo, B(r) « 1 — which represents 
Schwarzschild black hole, but it becomes de-Sitter 
black hole as B(r) « 1 — jr near the center (r ~ 0), 
so it is a regular space-time without singularity. Thus 


Hayward black hole is the simplest regular black hole. 
Some physical consequences of Hayward black hole 
have been discussed by several authors 40 - 142} . After 
that the Hayward metric was modifies [43} by choos¬ 
ing A(r) = f(r)B(r ), satisfying the conditions: it (i) 
preserves the Schwarzschild behaviour at large r, (ii) 
includes the 1-loop quantum corrections and (iii) al¬ 
lows for a finite time dilation between the center and 
infinity. So the modified Hayward black hole metric 
is given by 0 , 


ds 2 = —f(r)B(r)dt 2 + 


1 


B(r ) 


dr 2 + r 2 (d0 2 + sindd(j) 2 () 17) 


where 


B(r) = 1 - 


2 Mr 2 
r 3 + 2 Ml 2 


f(r) = 1 - 


a/3M 
ar 3 + j3M 


(18) 


with a, fj are positive constants. Now from the rela¬ 
tion A(r) = f(r)B(r), we may obtain 


A'{r) f'(r) B'(r) 
A(r) f{r) B(r) 


Also from the expressions of B(r) and f(r) (equation 
(HU)), we get 


B'(r) _ 2Mr{r 3 -4Ml 2 ) 

B(r) (r 3 + 2M 2 )[r 3 + 2M(1 2 - r 2 )] ^ 


f'(r) = _ 3 a 2 /3Mr 2 _ 

f(r) (ar 3 + (3M)[ar 3 + (1 — a)/3M] 


Since for outside the horizon, (i) B(r ) > 0 which 

implies r 3 > 2 M(r 2 — l 2 ) and (ii) f(r) > 0, so we get, 
1 

r > 3 w j^j 1 a > i g 0 f rom equation (l2ll) . 

we have f'(r ) > 0. 


If we assume that the fluid flow accretes upon mod¬ 
ified Hayward black hole, we can calculate the expres¬ 
sions of u 2 , V 2 and c 2 at the critical point r c . The 
expressions are given below (using equations (fT2l) . (fl3l) 
and flHJ) ): 


2 B'(r c ) (B'(r c ) 
2 V B(r c ) 



HMY 1 

f( r o) ) 


m) 


v 2 

y c 


/Bp) 
B'(r c ) V B(r c ) 


f'(r c ) V 1 (2 f(r c )\ 
f( r c) ) \r c f(r c )J 


, 2 = C A V c (V 2 + 1) 

u c ^f(r c ) 


(24) 


where B(r ), f(r) and their derivatives are given in 
(fl8l) . (EU1) and (T2TT) at the point r = r c . The physically 


(23) 
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acceptable solutions of the above equations may be 
obtained if u 2 > 0 and V 2 > 0 which leads to 

B'(r c ) > 0 and 0 < 4^4 < ~ (25) 

f\ r c) r c 

From above restrictions, we may get the bounds of r c 
and that is (a > 1): 


r 3 > Max 



/?(-4 + 5a + y/a(25a - 24))M 
4a 



For example, we assume a fluid flow obeys linear 
equation of state p — wp (w = constant) accretes 
upon modified Hayward black hole. Then we obtain 
c 2 s = w and V 2 = 0 and from (flOl) . we obtain u c = 0. 
From equations (1221) and (1241) , we see that the critical 
point occurs at the point r c = (4 Ml 2 )3. For general 
equation of state where w = w(t), we obtain c 2 ^ 
constant, V 2 ^ 0 and u 2 ^ 0. In this case, it is very 
difficult to obtain the critical point r c . 


IV. CHANGES OF BLACK HOLE MASS 
DURING ACCRETION AND EVAPORATION 

The rate of change of mass M of the black hole is 
computed by integrating the flux of the fluid over the 
2-dimensional surface of the black hole and given by 
M = — f TgdS where dS = y/^dddfy. Using equa¬ 
tion (O, we obtain the rate of change of mass of black 
hole as in the following form: 

M = 4tt CM 2 { Poo + p(p oo)) (27) 

The above result is also valid for any equation of state 
p = p(p)- So the rate of change of mass for the accret¬ 
ing fluid around the black hole will be 

M acc = 4 -kCM 2 {p + p) (28) 

We see that the rate of change of mass for the 
general spherically symmetric static black hole due 
to accretion of fluid flow becomes exactly similar 
rate in the case of a Schwarzschild black hole. From 
the expression (1281) it is to be noted that the rate of 
change of mass of any static spherically symmetric 
black hole is completely independent of A(r) and 
B{r). When some fluid accretes outside the black 
hole, the mass function M of the black hole is 
considered as a dynamical mass function and hence 
it should be a function of time also. So M is time 
dependent and the increasing or decreasing of the 
black hole mass M sensitively depends on the nature 
of the fluid which accretes upon the black hole. If 
p + p < 0 i.e., for phantom dark energy accretion, the 
mass of the black hole decreases but if p + p > 0 i.e., 
for quintessence dark energy accretion, the mass of 


the black hole increases. 

We may also assume that the black hole evaporates 
by Hawking radiation process. The rate of change of 
mass for the evaporation is given by 

Meva = (29) 

where D > 0 is a constant whose value depends on 
the model [44[. Now due to accretion of fluid flow and 
evaporation of mass of black hole, we get the rate of 
change of mass of black hole as 

M = M acc + M eva = 4 nCM 2 ( p + p )-JL ( 30 ) 

For accretion scenario, the changes of mass of black 
hole completely depends on the nature of the fluid 
accretes. But for evaporation process, the change of 
mass of black hole is independent of the nature of 
fluid, because this is internal process. In fact, when 
the accretion fluid is only the cosmological constant 
(p = —p), the mass of black hole for only accretion 
scenario is always same throughout time evolution. 
Only in accretion process, the mass of black hole 
increases for normal fluid and quintessence type dark 
energy fluid and decreases for phantom dark energy. 
But due to accretion as well as evaporation, M > 0 
for M4 > 4 *c? P +p) and M < 0 for M 4 < 4 „ c f p+p) for 
normal fluid and quintessence type dark energy, but 
for phantom energy, black hole mass always decreases 
(M < 0). Thus evaporation supports the decreasing 
of the mass of black hole with some restrictions of 
minimum values of mass of black hole. 


V. DISCUSSIONS AND CONCLUDING 
REMARKS 


First we have assumed the most general static 
spherically symmetric black hole metric. The ac¬ 
cretion of any general kind of fluid flow around the 
black hole have been investigated. For this general 
kind of static black hole, the critical point, velocity 
of sound, fluid 4 velocity have been calculated and 
shown that these value depend completely on the 
metric coefficients. Next, the accretion of fluid flow 
around the modified Hayward black hole have been 
analyzed and we then calculated the critical point, 
fluid 4 velocity and velocity of sound during accretion 
process. We can mention that for outside the horizon, 
(i) B(r) > 0 which implies r 3 > 2 M(r 2 — l 2 ) and (ii) 


f{r) > 0, so we get, r > 


P(a-l)M 

O' 
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with a > 1 


and also f'(r) > 0. For physical region of accretion 
have been found and the bounds of critical point 
have been generated and the bound of critical point 
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is r-3 > Max {4 Ml 2 , |. when 

the perfect fluid satisfying linear equation of state 
p = wp (w = constant) accretes upon modified Hay¬ 
ward black hole we have obtained c 2 = w, u c = 0 and 
V 2 = 0. In this accretion process, we have seen that 
the critical point occurs at the point r c = (4 Ml 2 )i. 
Also the nature of the dynamical mass of black 
hole during accretion of fluid flow and taking into 
consideration of Hawking radiation from black hole 
i.e., evaporation of black hole have been analyzed. 
Only in accretion process, the mass of black hole 
increases for normal fluid and quintessence type dark 
energy fluid and decreases for phantom dark energy. 


But due to accretion as well as evaporation, M > 0 
for M4 > 4 kC? p +p) and M < 0 for M 4 < 4 „ c f p+p) for 
normal fluid and quintessence type dark energy, but 
for phantom energy, black hole mass always decreases 
(M < 0). Thus evaporation supports the decreasing 
of the mass of black hole with some restrictions of 
minimum values of mass of black hole. 
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